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Abstract 

We provide a method for extracting information on the energy spectrum 
of solar neutrinos directly from the spectrum of scattered electrons. As an 
example, we apply it to the published Super-Kamiokande data. When combined 
with data from SNO on charged current interactions this method allows to 
derive separately the spectra of u e and of plus v T . 
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1 Introduction 



As a general rule, important information is contained in the spectrum of the radiation which 
is being used as a probe of a physical system. This holds for any radiation and thus also for 
neutrinos emitted from the Sun. Indeed, the reconstruction of the solar neutrino spectrum 
from experimental data is very interesting, as a way to establish neutrino properties and/or 
to study the stellar interior. 

Solar neutrino experiments provide however indirect information about the neutrino 
spectrum. As an example, from radiochemical experiments |TL |2[ ||] , which detect the decay 
of nuclei produced by neutrino capture over the target atoms, one can derive an average flux 
of electron neutrinos arriving onto earth, where the average is weighted with the capture 
cross section and it generally involves neutrinos from different branches of the fusion chain. 

Super-Kamiokande (SK) || is essentially sensitive to 8 B neutrinos only. It studies 
Electron Scattering (ES), 

v + e~ — ► v + e~ , (1) 

by measuring the energy distribution of electrons for total energies E e > 5.0 MeV. Electron 
scattering is being studied also by SNO, see |5), for kinetic energies T e > 6.75 MeV. 

ES measurements have given, so far, only indirect information on the distribution 
of the neutrino energy E u . In the usual approach, one starts with a flux (p a (E u ), where 
the index a corresponds to active neutrinos, and evaluates the scattered electron spectrum. 
If this is consistent (inconsistent) with experimental data then the input <p a is accepted 
(excluded) for describing the solar neutrino flux arriving onto Earth. 

The main purpose of this letter is to provide an alternative method for extracting 
the energy spectrum of solar neutrinos directly from data on the spectrum of scattered 
electrons. 

The basic idea is very simple. Electrons with total energy E e (which we assume much 
larger than m e ) are produced by neutrinos with E v > E e — m e /2. Above this threshold, the 
cross sections da a /dE e are practically independent of electron energy. Thus the difference 
of electron spectra centered at two neighbouring energies around E e gets contribution only 
from neutrinos with energy around E v = E e — m e /2. Consequently, the neutrino spectrum 
can be calculated from the derivative of the measured electron energy spectrum. 

First, we shall present this scheme by using simple analytical approximations. Next 
we shall confirm these estimates by means of a full numerical calculation. 

As an example, we shall apply this method to the published SK data see our 
results in fig. [5|. We recommend however that a full analysis be done by the experimental 
group, since a detailed knowledge of the detector and of the data is important for extracting 
optimal information. 

From ES experiments one derives information on a specific combination of the spec- 



2 



trum of u e , (p e , with the spectrum of plus u T , (f^ T : 

PES =Ve + @<Pij,t (2) 

where j3 is a suitable ratio of Vn (u T ) to v e cross sections. An important addition of SNO 
||], by means of the Charged Current (CC) data on deuterium, 

v e + d —> p + p + e~ , (3) 

is the possibility of determining the v e spectrum. We will show that, by combining ES and 
CC data, one can determine separately the spectrum of active neutrinos different from v e . 

We remind that the comparison of the (energy integrated) electron signal from ES 
and CC data has already allowed to disentangle the presence of a v^ T component in the 
electron scattering data |5|, g, 0]. By the method which we are proposing it becomes possible 
to determine the v„ T energy spectrum separately. This is not only interesting by itself, but 
it also offers a possibility to find a model-independent signature of sterile neutrinos in case 
they too are coming from the sun, as we shall comment at the end of the paper. 



2 From electron spectrum to the neutrino spec- 
trum 

a) A simplified case 

In the limit of infinite energy resolution and assuming full efficiency for electron detection, 
the electron energy spectrum from reaction (|l|) is given by: 



dN «r , , • 

dE P = Mt ldE » 



(4) 



where M is the number of target electrons, t is the measurement time, <p a are the neutrino 
fluxes and da a /dE e are the differential cross sections for active neutrinos (o = e,/i,r) [Q, [9|]. 

For energies much larger than m e , electrons with total energy E e are produced by 
neutrinos with energy larger than E m [ n = E e — m e /2. In addition, to a good approxima- 
tion in the range of measured electron energies, the differential cross sections da a /dE e are 
independent of the energy of the scattered electron. In other words, we can write: 

do a 



dE P . 



a a {E v ) 6{E U - E min ) , (5) 



where the factors a a are, moreover, weakly dependent on the neutrino energy E v . Therefore, 
by differentiating both sides of eq. (||) with respect to the electron energy we obtain: 

/ i \ d 2 N 

^) + ^ T ( E „) = -( i _j_. ( 6 ) 
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The factor (3 = a^/a e is weakly dependent on E u in the energy of interest to us and the 
r.h.s. is calculated at: 

E e = E u + m e /2 . (7) 

Eq. (^), which is our basic result, is the direct information on the neutrino energy 
spectrum which can be derived from the electron spectrum. We remark that the r.h.s is 
purely defined in terms of measurable quantities. 

Of course, there is no way at this stage to tell which are the separate contributions of 
v e and z^ T . This separation requires an additional information, such as it can be provided 
by the CC measurement of SNO, see below. 

On the other hand equation @ allows a direct test of oscillation models. Any pro- 
posed oscillation solution predicts a definite expression for the produced flux tp(E u ) and for 
the oscillation probabilities P ea {E v ). In terms of these, the l.h.s of eq. @ can be immediately 
calculated as <p [P ee +/3(P eiJj +P eT )] and compared with the observable quantity on the r.h.s. . 

b ) The effect of the finite energy resolution 
In practice, the situation is more complicated, since the experiment has a finite energy 
resolution. This implies that the observed electron energy e e is different from the true 
electron energy E e . The observed spectrum S(e e ) is related to the true energy spectrum 
dN/dE e by means of the following relation: 

S(e e ) = J dE e r(e e ,E e )^- (8) 

where the resolution function can be taken as a gaussian: 

By deriving both sides of eq. (|8|) with respect to the observed energy e e we obtain: 

™ = f dEe ?^IA™, (10) 

de e J de e dE e 

If we neglect the energy dependence of A , from @ we have dr/de e = —dr/dE e . Integration 
by parts gives: 

M _ ^,^N{E e ) 
de 

By using eq.@ we obtain: 

dS f 

— = -Mta e J dE v r(e e ,E u + m e /2) [(p R {E v ) + /3<p^(E v )] , (12) 

where we have taken advantage of the weak energy dependence of a e (E u ) to take it out 
from the sign of integration. 

Since experimental results are generally presented in terms of the ratios to the SSM 
prediction, it is convenient to write eq. ([l2] ) in a slightly different form. 



JdE e r(e e ,E e ) d -^l. (11) 
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-Nta e J dE v r(e e ,E v + m e /2) V SSM (E v ) . (13) 



We note that the electron spectrum predicted by Standard Solar Model (SSM) cal- 
culations, dS^ SSM ^ / de ei satisfies an equation similar to (|T2|): 

de e 

Also, for each energy we can normalize the neutrino fluxes to the SSM prediction, tp SSM , 
by introducing the quantities: 

<Pa(E; 



UE V ) = • (14) 



With these definitions, from eqs flllq ) and (13) we have: 

J dE v \f e (E v ) + PUr{E v )\ p(e e , E v ) = d / { ssM) e /dte ( 15 ) 

where the response function p(e e ,E u ) is given by: 

/ w s r(e e ,E u + m e /2) V SSM (E U ) 
p(e e ,E u ) = — r . (16) 



J dE v r(e e ,E u + m e /2) ip SSM \E V ) 



The function p(e e , E u ) essentially measures the energy resolution which can be attained 
for determining the neutrino spectrum (or, more precisely, its deviation from the SSM 
prediction). 

Eq. (|l5|) , which is our main result, is the extension of (||) for the case of finite energy 
resolution. It has a natural interpretation: due to the finite energy resolution, the derivative 
of the observed electron spectrum determines the neutrino spectrum, smeared over the 
energy resolution. 

The relation between neutrino and electron energy, previously given by eq.(0), can 
now be derived as follows. For a fixed electron energy e e , the l.h.s of eq.(|l5|) receives 
contribution from neutrino energies around E v such that dp/dE u = 0. Neglecting again the 
energy dependence of A this gives: 

X] lo SSM 

e e = E v + m e /2 - A 2 dh )f . (17) 



c) The general case 

We remind that equations ( |i~5D , (|l^) and (|17]) have been obtained by neglecting the electron 
energy dependence of da a /dE e and of A. If these approximations are released one obtains: 

dS/ de e 
dSiSSMT/de e 

where one has now two response functions: 

dr(e e ,E e ) da a (E u ,E e 



J dE v \U p e + Pfi] = JsTl) e /de 



V SSM (E U ) J dE e 



P., i < , . E„ ) — ^ . (19) 

^ a 
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The normalizations ^ a are given by: 

* = fiB. ^(B.) JdE c ^EA d ^S^A . (20 ) 
The factor f5 is expressed as: 

0(e e ) = ^ ■ (21) 



We can show that for SK eqs. (|15[|17|) are quite accurate. In fact, in the case of SK, 
the energy resolution can be described by the expression (||) with A given by 1 1 ] : 



A = 1.5 MeV^(E e - m e )/10MeV (22) 

By using this relation, the exact expression for the differential cross sections da a /dE e || 
and the SSM neutrino flux y^^Q, one can calculate numerically the response functions p e 
and for Super-Kamiokande. In fig. |l| we present our results as a function of neutrino 
energy, for selected values of e e . The functions p e and p M are bell shaped functions, with a 
full-width-half maximum of about 2.5 MeV. They are slightly narrower than the electron 
energy resolution r, due to the presence of (p SSM in their definition. 

One sees that, a part for the smallest electron energies, the exact response functions 
p e and are close to the approximate expression p given by eq. (|l6|) . For our purposes we 
can thus safely assume: 

p e (e e , E v ) = p^ee, E v ) = p(e e , E u ) . (23) 

In fig. U we show the relationship between neutrino and electron energies, calculated 
numerically by solving dp a /dE y = 0. One sees that the approximate solution given by eq. 
( |l7|) is quite accurate. 

Finally, in fig. || we present an estimate for /3, calculated according to eq. (^T|) , as 
a function of the electron energy. As the energy increases, it becomes approximatively 
constant, (3 ~ 0.15. 



3 Extraction of the v energy spectrum from SK 
data 

As an example, we outline here a numerical procedure for extracting physical information 
on f a (E u ) = ip a /ip SSM from the published SK data [f§ by means of eqs. (p^-17). 

*We remark that the response functions depend only on the shape of the neutrino spectrum. SK 
is essentially sensitive to 8 B neutrinos only, with a small contribution from hep neutrinos. We use 
the 8 B neutrino spectrum given in fl(i|| , and the ratio between 8 B and hep neutrino flux predicted 
by fl^l . Variations of the hep neutrino flux within the current phenomenological limits do not affect 
our results. 
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We remind that SK measures the electron spectrum S(e e ) for electron energies in the 
range 5-20 MeV, see fig. |3]. Data are grouped into 18 bins, each 0.5 MeV wide, covering 
the range 5-14 MeV, with an additional bin extending from 14 to 20 MeV. The signal in 
the i-th bin, Si, is associated with an error ASi, which we take (for the moment) as the 
sum in quadrature of the statistical and uncorrelated systematical errors given in @. Data 
correspond to 1258 days of exposure All quantities shown in fig. |3] are normalized to 
the predictions g^ SSM ^ corresponding to the Standard Solar Model (SSM) of [T^f] 

In order to evaluate the r.h.s. of eq.(|i~5|) and the associated error, it is convenient 
to divide the energy range probed by SK in a few intervals, each with a width comparable 
to the neutrino energy resolution, since the experiment is anyhow unsensitive to spectral 
deformations on smaller scales. On the other hand, by considering intervals larger than the 
original SK binning it is possible to decrease statistical fluctuations in the evaluation of the 
spectrum derivative. 

We have considered five intervals, with energies in MeV between (5-7), (7-9), (9- 
11), (11-13) and (13-20) respectively. Inside each interval S(e e ) and S^ SSM \e e ) can be 
approximated by straight lines and their slopes are determined by least squares fitting: 



dS _ (Se e )-(S)(e e ) 



de e (e 2 e ) - <6e> 2 



(24) 



and similarly for dS^ SSM ^ /de e . In the average ( ) each point is weighted according to the 
experimental error ASi: 

, ys _ Ej Xj/ASj 2 

For each interval we get 

dS/de e (Se e ) - (S)(e e ) 



dS (SSM) /dee (S(SSM) e J _ (S&SM)){ 
with an associated error: 

[(4) ~ (e e ) 2 ] 1/2 I r 



{S (SSM) ee) _ {S (SSM) ){ee) V ^ 1/AS 2 



(27) 



At this point we can calculate, by means of eqs. (15-17]) the deviations of the neutrino 
spectrum from the SSM predictions by using the SK data. 

Our results are presented in fig. |as a function of the neutrino energy. The quantity 
on the vertical axis represents f e + /3/ MT averaged with the response function p, i.e.: 

Fes(E„) = J dE' v \fe(K) + PMK)] p{te,K) (28) 

t The SSM of (H) predicts a 8 B neutrino fluxes $ B = 5.05 • 10 6 cm" 2 s _1 . Here and in the 
following, we use this value as a (convenient) normalization factor. We remark however that our 
results do not depend at all on solar models. 
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where e e (E u ) is given by eq.(|l7]) and (^). The horizontal bar corresponds to the neu- 
trino energy resolution calculated as the full width half maximum of the response function 
p(e e ,E u ). 

Concerning errors, the inner vertical bar takes into account statistical and energy 
uncorrelated systematical errors. The outer bar also accounts for energy correlated system- 
atical errors, as given in Q. Their effect has been evaluated by a simultaneous up and down 
shift of all the experimental points. 

We remark that from SK we have been able to derive only a specific combination of 
v e and Ufj, T fluxes, given by eq.p8|). For extracting the individual contributions of v e and 

ZV' i- e " 

F a (E u ) = J dE' u f a (E' u ) p(e e ,E' u ) (a = e, p,, r) (29) 
additional information are needed. 



4 The information from v e + d ^ p + p + e 

SNO has recently presented || the energy spectrum of electrons from reaction (||). We 
show here that these data, which provide a direct determination of the electron neutrino 
spectrum, can be combined with ES data for determining the spectrum of plus v T . 

The SNO results are shown in the lower panel of fig. | as a function of the observed 
electron kinetic energy T e = e e — m e . Data are grouped in 11 bins covering the energy range 
T e = 6.75 — 13 MeV. The signal in each bin, Cj, is normalized to the predictions, Cf , 



corresponding to the SSM of [12] and the error bars take into account only statistical errors. 

The measured electron spectrum C{e e ) from reaction (|3|), normalized to the SSM 
expectation C (e e ), is related to the the electron neutrino spectrum by: 



C 



J dE v f e (E u )p cc (e e ,E u ) (30) 



(jSSM 

where p cc is a suitable response function, defined in terms of the cross section for reaction 
(|3|), da cc /dE e [O], and of the detector resolution function rs^o(e e: E e )^. 

d<J C c{E u , E e ) 



^SSM 



(E v ) J dE e r S No(e e ,-E , e 



Pcc(e e , E v ) = J - ^ . (31) 

™ cc 

The normalization factor ^ cc is given by: 

* cc = J dE u <p SSM (E v ) J dE e r SNO (e e , E e ) da ^E e ) ^ 

The behaviour of the SNO response function is shown in fig.^ as a function of neutrino 
energy, for representative e e values. One sees that for each measured electron energy e e the 



* The SNO resolution function can be described by rel. (g) with A = (-0.4620 + 0.5470 *jE~ e 
0.008722 E e ) Mev [|. One can easily check that r SNO ^ r. 



s 



response function is peaked at a specific neutrino energy E v . The rates Ci/Cf can thus 
be interpreted as a determination of the deviations of the electron neutrino spectrum at E u 
averaged over the response function p cc , i.e.: 

F SNO (E u ) = J dE' v f e {E' u )p cc {e e ,E' u ) (33) 

where the numerical relation e e (E u ) is shown in fig. 

We remark that the width of the SNO response function, p cc , is mainly determined 
by the SNO resolution tsno- in fact, if one neglects the recoil of the two protons in 
the final state, there is a fixed relation between the neutrino and electron energies, i.e. 
da cc /dE e oc 8{E e — E v — nid + 2m p ). As a consequence one expects, roughly: 

p cc (e e , E v ) oc ip SSM (E u ) r SN o(e e , E v + m d - 2m p ). (34) 



We note that the structure of the previous relation is similar to that of rel. (16), 
which defines the response function p associated to the SK spectrum derivative. If one 
considers that SK and SNO energy resolutions are almost equal, this suggests that SK and 
SNO response functions can be equalized with a proper choice of the detection energies. 

The qualitative argument given above is clearly oversimplified. However, the possibil- 
ity to equalize the SK and SNO response functions can be checked numerically. The results 
are shown in fig. [7|. One sees that the response function p e and p cc , calculated numerically 



according to rel. ( |19| ) and (31) are almost equal, i.e. 

Pe{esK,E u ) = pcc (esno,E v ) (35) 
if one chooses the energies as follows: 

e SNO = 0.975 e SK - 2.50 MeV . (36) 



From the equalities ( p3| ) and (^5[), it follows that one can identify the SNO charged current 
spectrum with the v e contribution to the SK spectrum derivative, i.e.: 

F e {E v ) =F SNO (E u ) . (37) 

As a consequence, SK and SNO-CC data can be combined to determine: 

F„ T (E V ) = i (F ES (E U ) - F SNO (E u )) (38) 

and: 

F e (E v ) + F„ T (E V ) = i [Fes(E v ) - (1 - P)F SNO (.E u )] . (39) 

The results obtained are shown in fig. || where, for convenience, the SNO-CC signals have 
been grouped into larger intervals so as to reduce statistical fluctuations and to take advan- 



tage of the SK-SNO correspondence, eq. (|36j). 
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We remark that, if there are no sterile neutrinos, then the sum of the spectra of active 
neutrinos has the same shape as the 8 B spectrum in the laboratory, i.e. 

ip e {E v ) + ^ T (E U ) = k ip SSM {E v ) (40) 

where A; is a normalization constant, i.e. it does not depend on E v . Violations of this sum 
rule, would clearly imply sterile neutrinos. Such deviations are not shown in the data, which 
are well consistent with a constant, see the lowest panel of fig.pl 



5 Concluding remarks 

We have provided a method for extracting information on the neutrino spectrum from the 
spectrum of scattered electrons, summarized in eq. (15-17]). We summarize here a few 
points, concerning the possible applications and developments of our approach: 

1 ) As an example we have applied our approach to the published SK data, see fig. [5|. 
We suggest that the analysis is performed by the experimental group, since a detailed knowl- 
edge of the detector is important for extracting optimal information. 

2) We have shown that the information obtained by using our method can be directly 
combined with the information provided by SNO, so as to determine the spectrum of plus 
v T as well as the total spectrum of active neutrinos. Our results, obtained from eqs. ( [38]) 
and (|3^), are shown in fig. |j| 

3) We remark that the SK-SNO comparison potentially allows for a model indepen- 
dent signature for sterile neutrinos. If there are no sterile neutrinos, in fact, the sum of the 
spectra of active neutrinos has the same shape as the 8 B spectrum in the laboratory, i.e. 
ip e (E u ) + (furiEn) = k ip SSM (E u ), where k is a normalization constant. Violations of this 
sum rule, would clearly imply sterile neutrinos. Such deviations are not shown in the data. 
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Figure V. Response functions of Super-Kamiokande as a function of neutrino energy 
E u , for selected values of the observed electron energy e e . The exact response functions 
p e and pp, are calculated according to eq . (|T9|) , whereas p is the approximation given 
byeq.®. 
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Figure 2: The relation between neutrino and electron energies for SK and SNO de- 
tectors. The upper lines refer to SK. They are obtained from numerical calculation 
(solid) and from the approximate relation ([17|) (dot-dashed). The dashed line refers 
to SNO. 
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Figure 3: The function {3 from eq. (pl|) . 
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Figure 4: Upper panel: The electron energy spectrum measured by SK 0] normalized 
to the SSM prediction [12] . Lower panel: The electron energy spectrum measured by 
charged current reaction (Q) at SNO |J normalized to the SSM prediction |T2 . 
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Figure 5: Deviations of the neutrino spectrum from the SSM prediction, as a function 
of neutrino energy E u , from SK data. The horizontal bar is the neutrino energy res- 
olution. The inner vertical bar denotes the uncorrelated error, the outer bar includes 
the effect of energy correlated systematical errors. 
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Figure 6: The SNO charged current response function p cc , eq.([H]), as a function 
neutrino energy E u , for selected values of the observed electron energy e e . 
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Figure 7: a) The SK response function p e , eq. (|T9D, for the indicated electron energies 
e SK (dot-dashed lines), b) The SNO response function p cc , eq. flSTD , for the electron 
energies esNO given by eq. (^6|) (solid lines). 
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Figure 8: Deviations of the neutrino spectra from the SSM prediction as a function of 
the neutrino energy E v . The vertical bars take into account only statistical errors. In 
the upper panel, the full circles correspond to Fes, while the open circles correspond 

to Fsno 
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